Decoherence and dephasing errors caused by D.C. Stark effect in rapid ion transport 



O 

(N 

u 

< 

oo 

(N 

S 1 

Oh- 
•*-> 

cr 



> 
in 
\o 
o 
in 

cn 
o 



X 



Hoi-Kwan Lau* and Daniel F.V. James^ 
Department of Physics, University of Toronto, 60 Saint George Street, Toronto, M5S 1A7, Ontario, Canada 

(Dated: January 18, 2013) 

We investigate the error due to D.C. Stark effect for quantum information processing for trapped 
ion quantum computers using the scalable architecture proposed in J. Res. Natl. Inst. Stan. 103, 
259 (1998) and Nature 417, 709 (2002) 0. 

As the operation speed increases, dephasing and decoherence due to the D.C. Stark effect becomes 
prominent as a large electric field is applied for transporting ions rapidly. We estimate the relative 
significance of the decoherence and dephasing effects and find that the latter is dominant. We find 
that the minimum possible of dephasing is quadratic in the time of flight, and an inverse cubic in 
the operational time scale. From these relations, we obtain the operational speed-range at which the 
shifts caused by D.C. Stark effect, no matter follow which trajectory the ion is transported, are no 
longer negligible. Without phase correction, the maximum speed a qubit can be transferred across 
a 100 micron-long trap, without excessive error, in about 10 ns for 40 Ca+ ion and 50 ps for 9 Be + 
ion. In practice, the accumulated error is difficult to be tracked and calculated, our work gives an 
estimation to the range of speed limit imposed by D.C. Stark effect. 

PACS numbers: 03.67.Lx, 03.67.Pp, 37.10.Ty 



I. INTRODUCTION 



Since it was realized that quantum algorithms can 
speed up complicated computational tasks, such as fac- 
torizing a large integer, which cannot be efficiently per- 
formed by known classical algorithms, building a quan- 
tum computer (QC) has become one of the ambitious 
goals in modern physics [!, H[ . Among current proposals 
for physical implementations of a QC, in many ways the 
ion trap proposal of Cirac and Zoller [f| seems the most 
auspicious at the moment. Exploiting well-developed 
techniques from quantum optics and atomic physics, en- 
tanglement of up to 14 ions [f| 0], high fidelity gates and 
readout 0, |9| , and long coherence time quantum memory 
for more than 10 s [lfjj] have been recently demonstrated. 
Simple quantum algorithms including Deutsch-Jozsa al- 
gorithm [ll[ and quantum teleportation [l2|, [H| , as well 
as verification of the Bell inequality [3] have been suc- 
cessfully realized on ion-trap QC. 

Nevertheless, we are still far from having a quantum 
computer with the computing power higher than (or 
even comparable to) classical computers. Apart from the 
small number of qubits entangled, when comparing cur- 
rent capabilities with the resources needed for a practical 
quantum processor, the speed of quantum operations is 
another issue that limits the clock rate of an ion-trap 
QC. Consider a QC is built to run the Shor algorithm. 
It would be of great technical interest if it could break 
a RSA classical cryptography code in a few days, but it 
would obviously be even more appealing to break the 
code, for example, in a few hours. Then each quan- 
tum logic gate has to be performed at the time scale of 
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lis [T|| 0|. When employing fault-tolerant techniques, 
each quantum logic gate consists of several concatenated 
rounds of physical quantum operations on ion qubits, 
including transportation, cooling, and laser interaction 
[2, [3]. Assessing the true time specifications for a QC in 
a complicated problem dependent on paradigm choices 
such as circuit-based versus measurement-based QC, er- 
ror correction code used, and the implementation of the 
algorithm to be performed. But we can assert quantum 
computation is more useful and promising if the speed of 
each physical quantum operation is of nanosecond scale 
17]. 

In this paper, we investigate the effect of D.C. Stark 
shifts during ion-trap QC operation according to the scal- 
able architecture proposed by Wineland et al. [l| and 
further elaborated by Kielpinski, Monroe and Wineland 
in Q • The QC consists of large number of interconnected 
ion traps. Ions are stored in memory region, which is far 
from interaction region where entanglement operations 
take place. During quantum computing, ions are moved 
from the storage region to the interaction region, and 
transported back to storage region after operations. The 
transport of ions can be done by changing electric po- 
tential of each trap, so that an effective non-equilibrium 
electric field is induced. 

Fast transport of the ion requires a large electric field 
which will result in a variety of potentially detrimental 
effect. One problem is that a large electric field is less 
stable, the field fluctuations will heat up the ion to mo- 
tional excited states [ijj]. Cycles of sympathetic cooling 
are required to bring the ion back to its motional ground 
state for precise logic operations, but the operation is 
time-consuming and thus not preferable in high speed 
quantum computation. The motional heating effect is an- 
ticipated to be reduced by improving experimental tech- 
niques, such as using surface traps and coating the elec- 
trodes of the trap [l8[, or transporting the ions under 
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traj ectories with minimal vibrational quanta excitation 
19]. Besides, some proposals for entanglement gates re- 
main effective even though the ion has small motional 
excitation (20l,[2l|. 

More seriously, a large electric field will induce a 
D.C. Stark effect onto the internal electronic states of 
ions. Due to the detuning of energy levels and mix- 
ing of eigenstates, encoded quantum information will be 
phase-shifted and excited out of the computational basis. 
This issue was first discussed more than 10 years ago by 
Wineland et al. (see ref. p.310); its significance was 
neglected at that time due to the low operational speed 
they envisioned. As technology advances this effect will 
become important, and a reappraisal is warranted. The 
aim of this article is to find out the relation between the 
total effects by D.C. Stark shift, the size of trap, trajec- 
tory and total time of the flight of ion qubits. From the 
relation, we can define the 'threshold speed', the speed 
of ion-shuttling above which the influence of D.C. Stark 
effect becomes significant. 

Our paper is organized as follows. In section[Hl we cal- 
culate the wave function of the ion during the shuttling 
process. A quantum mechanical treatment is employed, 
and we find the ion remains in coherent state during the 
trajectory. Therefore the classical acceleration of the ion 
can be obtained, from which the net electric field expe- 
rienced by the ion is determined. The total dephasing 
caused by D.C. Stark shifts is calculated in section Hill 
In section IIV[ we minimize the dephasing effect with re- 
spect to the ion-transport trajectory; a threshold speed is 
deduced from the relation between minimum phase shift 
and the time of flight. In the next section [V] we calcu- 
late the threshold speed for various choices of ion qubits. 
Significance of the decoherence effect is discussed in IVII 
and we summarize our results in section IVIII with some 
discussion. For simplicity, we do not consider the junc- 
tions where ions are transferred between different traps; 
the ion is assumed to move along a straight trap only. 



II. MOTION OF ION 

Let us consider an ion is confined in a trap. Its de- 
grees of motion in the y and z directions are restricted 
by a strong radio frequency electric field, so the ion is 
constrained to move only along the x direction. In cur- 
rent experiments of ion trap QC, the potential of the r.f. 
trapping field is of the order of 100 V, oscillating at radio 
frequency at about 100 MHz [18]. In spite of the strength 
of the r.f. field, the wave function of the ion is concen- 
trated and peaked at the saddle point of the potential 
in the y — z plane, on which the electric field vanishes. 
Therefore the r.f. field only induces a quadrupolar Stark 
effect on the ion, which is a higher order effect than the 
dipole Stark effect during the shuttling and will be ne- 
glected in our discussion. For the same reason, the D.C. 
Stark effect caused by micromotion of the ion in the y — z 
plane is also neglected. 



In the x direction, the ion is weakly trapped by an elec- 
trostatic field. Close to a minimum of the axial electro- 
static potential, the Hamiltonian of the system can be ap- 
proximated by a harmonic oscillator (Here we only con- 
sider harmonic potentials, more general trapping poten- 
tial and detailed dynamics of trapped ion are discussed 
m During shuttling, the strength of the electro- 

static field is changed so that the potential well is dis- 
placed. Carefully tuning the strength of the electrostatic 
field can ensure the strength of harmonic well is constant 
throughout the process. Here we assume the electric field 
is well controlled so the motional heating caused by field 
fluctuation can be neglected, this assumption will be dis- 
cussed later. The time-dependent Hamiltonian is given 
by the formula 
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(1) 



where p is the momentum operator of the ion in x direc- 
tion; uj is the harmonic frequency of the potential well 
formed by electrostatic fields; s(t) is the position of the 
center of the potential well, which is time dependent . 

One can solve for the ion wavefunction ]*&(£)) governed 
by the Hamiltonian Eq. {jXJ) in a few steps. First, the ion 
should be oscillating around the center of the potential 
well, so we consider a wavefunction translated by a dis- 
tance s(t). This can be done by applying a displacement 
operator D on the state. We define a new state |x(i)) by 
the formula 



\x(t)) = D 



2H 



s(t) 



!*(*)> > 



where [22| 

D[v] = exp(wa^ — v*a) . 
The Schrodinger equation becomes 

ifid t \x(t)) = (H + s(t)P) \x(t)) 



(2) 



(3) 



(4) 



where Hq is the untranslated harmonic oscillator Hamil- 
tonian, i.e. Eq. (JlJ with s(t) = 0. 

In the interaction picture, i.e. exp (^—iHot^ \x(t)) = 

\x(t)), the wavefunction satisfies the following equation 

d t \m) = (ma-u*(t)tf)\m) > (5) 

where a and are lowering and raising operator of Hq; 
u(t) is defined as follows: 



u(t) 



mui 



2h J Q 



s(t') exp(iut')dt' 



(6) 



The solution of Eq. ([5]) is straightforward [22| : 

\m) = D[u(t)}e*W\ X (0)) , (7) 
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where <f>(t) is a time-dependent overall phase of the ion 
given by 



[u(t')u*(t') - u*{t')u(t')]dt' 



(8) 



Combining these results, the wavefunction of the ion is 
therefore given by the formula 



D 

xD 



exp (^—iHat/hj 



u(t)-\/—s(0) 



3* w |*(0)) . (9) 



If we assume the initial state of ion is the motional 
ground state of the harmonic oscillator centered at x — 0, 
i.e. 1^(0)) = |0), then the motional state of the ion 
remains in a coherent state throughout the whole process, 
i.e. = \a(t)), with the amplitude a(t) given by the 

formula 



«(*) = \/ 2ft (*(*)- ' 



s(ti)e dt\ 



The expectation value of position q{t) is hence 



q(t) = s(t)- s(ti) cos [w(t - ti)] dti 
Jo 



(10) 



(11) 



Because of the classical nature of coherent states, the 
net electric field, £ = ^e^ where e x is the unit vector 
along x direction, experienced by the ion is related to 
the acceleration of the expectation value of position by 
the Newton's third law, viz., 



-m 



e 

mcu 2 



s(ti) cos[u;(ii — t)]dti 



(12) 



III. PHASE SHIFT DUE TO D.C. STARK 
EFFECT 

In this section, we consider the total phase shifted dur- 
ing the shuttling process. The effect of the applied elec- 
tric field on the internal structure of the ion will be de- 
scribed by an additional Hamiltonian -ffstark = — d ■ £, 
where d is the dipole operator. Since the fastest shut- 
tling time achievable in foreseeable future would not be 
faster than 10 ps, which is about four order slower than 
the reciprocal of the atomic transition frequency, we can 
obtain the D.C. Stark shift energy at time t by the time 
independent perturbation theory [23j j . viz., 



\ei{t){m\x\n)^ 



frlJri 



where electric field points to z direction only; tkj nm = 
hio n — hu m is the energy difference between energy eigen- 
states \n) and |m) of the unperturbed Hamiltonian. 

Suppose two unperturbed energy eigenstates |/) and 
|z) are chosen as the computational basis of the ion qubit. 
We work in the interaction picture. The ion is initially 
encoded with some quantum information as the state 
a\i) + fi\f). Because the D.C. Stark energy of each com- 
putational basis states are different, the qubit state gains 
an extra phase after shuttling, 

a\i) + p\f) ^ a\i) + pe^\f) , 

up to an unimportant global phase. Substituting from 
Eq. (fT2")l and flT3"l) . the extra phase factor (j> is given by 



where 



C[q(t)} = / q(t) 2 dt 
Jo 



(14) 
(15) 



All terms except C depend on atomic structure of the ion 
only, so a good choice of ion qubit would give a small 
magnitude of these terms; £ is a functional of the trajec- 
tory of the ion q(t) which is the same for any choice of 
ions and any computational basis states being used. 



IV. MINIMUM POSSIBLE PHASE SHIFT 

From Eq. (IT4")) and (TT51) . we see that the total phase 
shift due to the D.C. stark effect is determined by the tra- 
jectory. Among all the trajectories the ion travel, there is 
one, the optimal trajectory q$, which produces the min- 
imum phase shift, (f> m in(L, T), depending on the length 
L and time T of the flight. When the ion is shuttled 
in a particular time scale at which m ; n becomes signif- 
icant, the phase shift has to be corrected by some mea- 
sures, either unitaries or quantum error-correcton-codes 
(QECC), regardless the trajectory's form and whether it 
is predictable or not. This thus places a threshold on 
ion-shuttling speed. 

To find the minimum possible phase m i n , we have to 
find an ion trajectory q(t), for which £ is a minimum. 
We solve this problem using calculus of variation (see, 
e.g. [13] )• By defining 



lit) = «(*), 



(16) 



(13) 



minimizing £ is equivalent to minimizing £, defined as 

C[<7,?,7.7.*]= / {j 2 (t) + »(t)[q(t)- 7 (t)]}dt , (17) 
Jo 

where fi(t) is the Lagrange multiplier corresponding to 
our assumption Eq. (|16p . It is beneficial to minimize £ 
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because it is a functional of q, 7, and their first order time 
derivative only (i.e. no high order derivatives). Then we 
can easily write down the Euler equation with respect to 
q and 7 



fi(t) - 2 7 (t) = 

m = 



(18) 
(19) 



Incorporating these equations and the constraint 
Eq. (Til?]) , we find go obeys 



%'(*)= 



(20) 



We assume the ion is displaced by distance L after 
the process, the initial and final position can be set as 
q(0) = and q(T) — L. Additionally, we require the ion 
remains in the motional ground state before and after the 
shuttling, so the initial and final velocity both vanish, i.e. 
q(0) — q(T) — 0. Putting these initial conditions into 
Eq. (|2"0")) , we find the optimal trajectory is 



(21) 



which gives C(?o) 
shift is 

12m 2 L 2 



\2L 2 /T 3 . Hence the minimum phase 



HT 3 



\(m\x\z)\ 2 



\(m\x\f)f 



fin 



(22) 

We note that the ion does not always stay at the equi- 
librium position of the potential well s(t). As expected 
by classical intuition, the ion keeps sloshing as the well 
is displaced. Transporting the ion in a desired trajectory 
has to be achieved by tuning the trapping electric field 
carefully to balance this sloshing. The position of the 
well can be obtained by combining Eq. (|11[) with its sec- 
ond derivative. The displacement of the well sq (t) which 
produces the optimal trajectory of ion is given by the 
following expression: 



*o(t) = qo(t) + 



qo(t) 



= L -2 



'T 3 



t 



6 



2 T 3 (uTf 



(23) 



During the shuttling, i.e. < t < T, and the ion should 
stay at the equilibrium point of the well before and after 
the travel, i.e. s(t < 0) = and s(t > T) = L. Both 
qo(t) and So(t) are shown schematically in Fig. [1] 

We can see the position of the potential well jumps 
sharply at t = and t — T . This means that the trapping 
electric fields are sharply changed at the beginning and 
at the end. However, the sudden variation of electric 
field does not cause a significant D.C. Stark effect. It 
is because the energy perturbation only depends on the 
strength of electric field but not its time derivative within 
the dipole approximation. 
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FIG. 1: Time dependence of optimal trajectory of ion qo(t), 
and optimal trajectory of potential well so(t) for different 
trapping frequencies. 



We have found the optimal trajectory which the phase 
shift is minimum, but the motion of ion cannot be ex- 
actly controlled or predicted in practice, and trajectories 
other than qo(t) may be employed due to experimental 
convenience. To verify the appropriateness of our thresh- 
old speed deduced from 4> m im we have to see if the phase 
shift would vastly increase if the ion takes other trajecto- 
ries. Firstly we consider a more experimental realizable 
trajectory which the electric field is tuned gradually, i.e. 
no sharp slips of potential well. We try the trajectory 
described by 



q{t) = L{6t 5 /T 5 - 15< 4 T 4 + 10t 3 /T 3 ) 



(24) 



where the position of the electric potential well varies 
continuously as time, even at the beginning and the end. 
We find for this trajectory that C = 17.1L 2 /T 3 , which 
is larger than C(?o) but remains in the same order of 
magnitude. 

Secondly, we consider the trajectory of the experimen- 
tal setting of Rowe et al. [l4[ , where the location of the 
trapping potential sr varies as 



2 / \ 



suit) = Ls'm ( — j 



(25) 



and the frequency of the potential well is lor — 2-7T x 
2.9 MHz. The trajectory of the ion qu can be obtained 
from Eq. (JTTJ) , and the value of trajectory dependent term 
C in the expression of phase shift Eq. (IT4")) is C{qr) = 
24.3L 2 /T 3 , therefore the dephasing is double of c/> m i n - 
For these two examples, the phase shifts are at the same 
order as <fi m i n , therefore the threshold speed derived from 
4> m in is a useful benchmark. 
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V. THRESHOLD SPEED OF TRANSPORTING 
ION QUBITS 

In this section, we investigate the threshold speed of 
common ion qubit candidtates. We consider what if the 
phase shift due to D.C. Stark effect is neglected, i.e. no 
unitary correction is conducted. It will be a problem if 
its magnitude is to large to be corrected by QECC of 
the system. The threshold error rate of the best known 
QECC now is about 1% (e.g. [25)]). To estimate the 
order of magnitude of the minimum shuttling time, we 
let the upper limit of the optimal phase in Eq. (|2"2"j) as 
7r/100. We assume the length of the trap being travelled 
is L = IPO/mi, which matches setting of current experi- 
mental setup [Tj]. Two popular proposals of ion qubits 
are investi gate d. One proposal uses electronic states of 
40 Ca + ion [12] as qubits levels. The computational ba- 
sis states are \i) = \Si/ 2 , Mj = —1/2) ground state and 
the I/) = \D 5 / 2 ,Mj = —1/2) metastable state, which 
the T\ time is about 1 s. Anoth er p roposal uses hy- 
perfine states of 9 Be + ion [13, [TH, [I3 |. The computa- 
tional states we consider are \i) = \F = 2, to = 0) and 
|/) =\F = 1, Mp = 1). The T-2 time is greater than 10 
s when the beryllium ion is subjected to magnetic field 
with strength 0.01194 T 10], while the Titime for these 
hyperfine states is longer than centuries {26J . 

Now we calculate the term inside the bracket in 
Eq. (|2"2")l . Values of the matrix elements can be calcu- 
lated from tabulated parameters found in ref. [27|. We 
calculate the term belonging to \i) while the one of |/) 
can be calculated by the same method. Let the state 
\i) belongs to energy level k, it can be expressed as a 
superposition of states with definite magnetic quantum 
number |r), 



(26) 



where g k is the degeneracy of the level k. The last re- 
lation is derived from the orthonormality condition of 
Clebsch-Gordan coefficients. 

For the calcium qubit, only the states \m) = 
\P 3 /2,Mj = -1/2) and \P 1/2 ,Mj = -1/2) can yield 
non-zero dipole matrix elements. By using the formal- 
ism mentioned above, we find the minimum phase is 



iCa 



9.86 x 10" 



IT 3 } 



(29) 



where the square bracket denotes the value of quantities 
in S.I. unit. Both computational states of beryllium qubit 
are hyperfine states of the ground electronic state, so 
|m) are states on P 1 / 2 and P 3 / 2 levels only. We find the 
minimum phase is 



mm 



2.6 x 10 



[T 3 ] 



(30) 



For L = 100 fim and </> m ; 
of flight for ion qubits are 



< 



7r/100, the minimum time 



T, 



Ca 



> 14.6 ns 



T^L > 0.044 ns 



(31) 
(32) 



T®f n is 3 order smaller than because the Zeeman 
energy between two qubit states, Ez, is small compared 
to energy difference between atomic energy levels, Ea- If 
no magnetic filed is applied and the tiny hyperfine split- 
ting is neglected, the value of the bracket in Eq. (|2"2"j) van- 
ishes since both qubit states are atomic ground states. 
The first non-vanishing term would be the first order 
term of E z /Ea ~ 10~ 6 , and so the phase change of 
beryllium qubit is much smaller than calcium qubit. 



We pick the intermediate states |m) as energy eigenstates 
with definite magnetic quantum number, so the summa- 
tion of | to) involves summing of states in a particular 
energy level I and then sum over all energy levels. The 
transition energy hu>i m is the energy difference between 
levels k and I and given explicitly in 27] . For the ma- 
trix elements of dipole operator, we extract the reduced 
matrix elements. 



m\x\%) 



IWI^IIOPEEl^Ha 



(27) 



where Q 1 is the rank 1 irreducible tensor operator, 
| (k\ \Q 1 1 \l) | 2 is the reduced density matrix between energy 
levels I, C mp is the ClebschGordan coefficients between 
| to) and \p) [28|. The reduced density matrix can be cal- 
culated from the line strength Sik between energy levels 
I and k provided in (27| . By definition, 

S lk = g k J2(p\d\™) ■ (m\d\p)=e 2 g k \(k\\Q 1 \\l)\ 2 , (28) 



VI. DECOHERENCE 

Apart from dephasing effect, the electric field also ex- 
cites the electron into states outside the computational 
basis states; neglecting this effect will cause decoherence 
error on quantum information. We start by writing the 
D.C. Stark effect Hamiltonian in the interaction picture 



Vr(t) = -V rn q(ty urnt \r){n\ , 
e 



(33) 



where 



777 777 — » 

(r\H stalk \n) = —V rn q(t) = (r\d ■ x\n)q{t) . (34) 

e e 

Using time dependent perturbation theory, the first three 
terms of the Dyson series of the propagator Ui(t) in the 
interaction picture are as follows: 
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Ui(t) « J- —Vm f q(t'y" rnt 'dt\r)(n\ - ^^5>-'^'« /' 9(*")9(*')e i(w "'" t " +a, -' t ' ) ^"*»<n| . (35) 



/o Jo 



The degree of decoherence is determined by the ampli- 
tude of non-computational states at t — T . The two most 
significant terms are: the first order term with u r ^ Li n , 
because the states in the same energy level have the same 
parity and so V nn — 0; and the second order term with 
u) r = uj n such that U) rr / -\-LU r ' n = 0, otherwise there will be 
a fast oscillating term in the integral and the contribution 
of the term is reduced by an order of perturbation. 

Although the time integrals cannot be solved without 
the exact form of q(t), we can compare the significance of 
both terms by estimating their magnitude. After simple 
integration by parts, the integral in the first order term 
will be 



<i{ty< 



l dt = (g(T)e l """ J -<?(0) 

UJ rn 



/ q^y^dt 

Urn Jo 



(36) 



Since the typical order of atomic transition frequency ui rn 
is 10 15 s _1 and the timescale of ion transport is about 
10 -9 s, the last term on the right hand side can be ne- 
glected as it is 10 -(> smaller than the one in the left hand 
side [23|]. As we are only interested in the order of mag- 
nitude, we consider the bracket in the right hand side has 
roughly the same magnitude as L/T 2 . 

For the second order term, we can do a trick to estimate 
its value. Consider the contribution for \r) = \n). This 
corresponds to the phase shift that we have calculated in 
previous sections. Therefore we know the time integral 
can be approximated by i((q)/uj r ' n , which the value is 
at the order of 10L 2 /T 3 . We consider V r ' n of low lying 
energy states dominate, and the matrix elements of the 
states in the same energy levels are in the same order. 
The ratio of the first and second order term is estimated 
by 



heT 



first order 



second order 10|K-J£ 



(37) 



For a qubit moving in a 100 /im trap within 10 ns, the 
above ratio is about 10 -4 , which means that the second 
order term should dominate the decoherence effect. 

The magnitude of the decoherence error is defined by 
the fidelity of the quantum state after passing through an 
error channel. Consider the D.C. Stark effect has excited 
the electron to non-computational basis state |6j), 



I*) 



(38) 



so the decoherence error is given by J^i \ e i\ 2 - For a cal- 
cium qubit, the amplitude of the second order term in 



Eq. ([33]) is on the same order as the phase shift (j>, so 
the magnitude of decoherence error is about (j> 2 . If the 
calcium ion is transported in a speed such that the de- 
phasing error is 10 -2 , the decoherence error will become 
10~ 4 — 10~ 3 . Therefore the threshold speed of trapped 
calcium ion QC is determined by the dephasing effect, 
i.e. Eq. (|3l]). 

The second order decoherence term of beryllium qubit 
vanishes if hyperfine splitting is neglected. It is because 
both computational states consist of one orbital state, 
\l = 0,1711 =0), only, i.e. 

lo n . (I .3 -1. .11. 1 .3 1. .1 -1. 



\1,1)f = 




(39) 



-I- -) r l- ^ 

2 ! 2' 2 ;/| 2' 2 7 



(40) 

where the first and second number of the states denote 
their angular momentum quantum number and magnetic 
quantum number; the subscripts specify the type of an- 
gular momentum. Since the orbital state has no degen- 
eracy, D.C. Stark effect can only produce a phase shift. 
Without hyperfine shift, the phase shift on each compo- 
nent of the computational basis is the same, so there is 
no decoherence. 

Just as the dephasing effect, the second order deco- 
herence term of beryllium qubit is non- vanishing only if 
hyperfine splitting is included, but its magnitude is re- 
duced by a factor of E Z /E A « 10" 6 . Then according 
to the previous analysis, the first order decoherence term 
may become dominant. However, the first order deco- 
herence term can be suppressed by tuning the trajectory 
of the ion qubit. Recall we have assumed the bracket 
in Eq. (|36|) takes the value L/T 2 , which determines the 
order of magnitude of the decoherence. However, if we 
choose a trajectory that the initial and final acceleration 
of the ion qubit are both zero, then the bracket vanishes. 
Hence the first order decoherence effect will be reduced 
by a factor of (u; rn T) _1 , which is a part in 10 6 for T = 1 
ns. Although the optimal trajectory qo(t) does not satisfy 
this criteria, as shown in Fig. [TJ we can slightly modify 
qo(t). An example to doing this is to increase the ac- 
celeration from zero to qo (0) initially and decrease from 
qo{T) to at the end of the flight, both are conducted 
slowly comparing to u rn but fast comparing to T, and 
the rest of the time follows qo(t). We find such modifi- 
cation does not alter the phase shift much, with only a 
factor of order 1, but the first order decoherence factor 
is greatly suppressed. 



So the second order decoherence term of beryllium 
qubit remains dominant. A similar situation pertains 
to the calcium qubit. Although the amplitude of non- 
computational states is on the same order as the dephas- 
ing of computational states, the decoherence error is the 
square of the amplitude. Therefore, we can state with 
some evidence that threshold speed of beryllium qubit is 
determined by the dephasing error, which is smaller than 
the order of 100 ps as given in Eq. (1321) . 

VII. SUMMARY 

In this paper, we have estimated the magnitude of D.C. 
Stark shifts on quantum information during ion shuttling. 
A larger electric field should be employed to transport ion 
qubits faster for rapid QC operation, but it also induces 
more serious D.C. Stark effect on the computational ba- 
sis states, which shifts the relative phase between qubit 
states and excits the electron out of the computational 
basis. Magnitude of these effects are characterized by the 
trajectory of ion transportation, so the threshold speed 
can be estmiatcd by minimizing the trajectory. We find 
that dephasing is the most important source of error, 
its magnitude scales as square of the size of traps and in- 
versely cubic as the operational timescale. For a trapped- 
ion QC with characteristics length about 100 /xm, with- 
out correction the shortest time for each shuttling T m j n 
is about 10 ns if 40 Ca + ion qubit is employed. On the 
other hand, using hyperfine states of 9 Be + ion as qubit 
states will be more error- resistant. We find that the to- 
tal phase shift due to D.C. Stark effect is 6 order smaller 
than the implementation using calcium ion. Beryllium 
ion can be transported as fast as 50 ps before dephasing 
and decoherence effects become significant. 

In principle, the average Stark shift might be assessed 
by Ramsey interferometry, and then corrected by unitary 
transformations. In practice, however, the ion is under- 
going a complex trajectory involving acceleration and de- 



celeration, moving in straight lines, turning around bends 
and through junctions, disengaging the individual ion 
from the storage register and the logic trap; all of these 
effects will be too complicated to track and calculate the 
Stark shift accurately. The uncertainties will turn the 
shifts into dephasing and decoherence errors. Error cor- 
rection will perforce be needed, and the requirements of 
fault-tolerant quantum computing (in particular, ensur- 
ing the error be below some threshold) will place a speed 
limit on the operation of the QC. Magnitudes of the er- 
rors depend on the setup of the ion-trap system. In any 
case, our result is a useful reference to the speed limit. 
As a rough illustration, suppose the overall Stark shifts 
can be evaluated, with accurately tracking the trajec- 
tory, well-controlled electric field, and other very precise 
experimental techniques, up to 90% accuracy for a par- 
ticular ion-trap QC. The speed limit of this setup, due to 
the 10% uncertainty, can be evaluated by Eq. and 
(|5U| as roughly two times of the threshold speed. 

A possible way to reduce this lower bound is to reduce 
the size of ion trap. But this method is inefficient be- 
cause T m j n only scales at L 2 / 3 , and the reduction of size 
by half already require tremendous advance in engineer- 
ing techniques. Another method is to encode quantum 
information onto the decoherence free subspace of two 
calcium ions [HH, e.g. |0) -> \hf 2 ), |1) -> As 
long as the traveled trajectory by both ions are approxi- 
mately the same, the total dephasing effect on the logical 
qubit can be reduced by several orders |2j . On the other 
hand, there are alternative scalable ion-trap QC architec- 
ture which requires much fewer times (the measurement- 
based ion-trap QC, c.f. Qjlllil), or even no ion-shuttling 
(the ion-photon network model, c.f. [3ll - l33| ); so they are 
not seriously affected by the D.C. Stark effect. 
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